The probability of double-strand breaks in giant DNA decreases markedly as the 

DNA concentration increases. 
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DNA double-strand breaks (DSBs) represent a 
serious source of damage for all living things and 
thus there have been many quantitative studies of 
DSBs both in vivo and in vitro Despite this 

fact, the processes that lead to their production 
have not yet been clearly understood, and there 
is no established theory that can account for the 
statistics of their production, in particular, the 
number of DSBs per base pair per unit Gy, here 
denoted by Pi, which is the most important pa- 
rameter for evaluating the degree of risk posed 
by DSBs. Here, using the single-molecule obser- 
vation method with giant DNA molecules (166 
kbp), we evaluate the number of DSBs caused by 
7-ray irradiation. We find that Pi is nearly in- 
versely proportional to the DNA concentration 
above a certain threshold DNA concentration. A 
simple model that accounts for the marked de- 
crease of Pi shows that it is necessary to con- 
sider the characteristics of giant DNA molecules 
as semiflexible polymers to interpret the intrinsic 
mechanism of DSBs. 

As for DSBs caused by 7-rays in vitro experiments, 
its reported values scatter between 7.1 x 10 -9 Q and 
2.1 x 1CT 6 [H (expressed as the number of DSBs per base 
pair per unit Gy) . The scattered data may be attributed 
to significantly varying DNA concentrations and the sizes 
of DNA. Although gel electrophoresis has been the most 
widely used method to enumerate DSBs, it requires a 
high DNA concentration (more than several tens of fjM. in 
base pair units) [§-0] ■ Contrastingly, the single- molecule 
observation method employed in the present study is ap- 
plicable to systems with DNA concentrations as small as 
three orders of magnitude smaller than those used in gel 
electrophoresis [f|. In this article, we report the marked 
dependence of DSBs on the DNA concentration over a 
wide range of concentrations, indicating the significance 
for genome-sized DNA. 
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FIG. 1. Representative fluorescence images of DSBs 
for T4 DNA. T4 phage 166kbp DNA in a 20 mM phosphate 
buffer solvent was irradiated by 7-rays from a Co 60 source. 
The fluorescent dye YoYo-1 was used for visualization of the 
DNA chains. DNA molecules in an aqueous solution were 
elongated through the effect of a weak shear applied by a 
glass slip and fixed on a poly-L-lysine coated glass, a, Gy 
(natural length), b, 50 Gy. c, 100 Gy. 



Figure [T] exemplifies the fluorescence microscopic im- 
ages that we obtained of the DNA in solution on the glass 
substrate. We obtained the average length of molecules 
irradiated with a dose /, L(I), for a range of / and 
thereby evaluated the number of DSBs per base pair, 
denoted by Po (see the Methods section). 
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FIG. 2. Irradiation dose dependence of the number of 
DSBs per base pair. The number of DSBs per base pair, 
Po, as a function of the irradiation dose from the Co 60 source, 
/. The data points for four different concentrations of samples 
are plotted. 
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FIG. 3. DNA base-pair concentration dependence of 
the number of DSBs per base pair per unit Gy. Log- 
log plot of the number of DSBs per base pair per Gy, Pi, as a 
function of the DNA base-pair concentration, Cbp. It is seen 
that Pi is roughly constant for small values of Cb P (compare 
with the solid line of zero slope) and inversely proportional 
to C'bp for large values. The blue dashed line represents a 
least-squares fit to the data using equation ©. The best fit 
is obtained with N Tad q/V tot = 1.31 pM (100 kbp ■ Gy)" 1 and 
V/N hp = 1.38 (a*M)-\ 



Figure [2] displays the dependence of Po on /. In the 
range of DNA concentrations from 0.015 to 300 /j,M in 
base pair units, P for all samples showed linear depen- 
dences on I. Thus, Pi for a given Cb p , the DNA base-pair 
concentration, is uniquely determined as the slope of the 
corresponding line representing Po(-0- 

Figure [3] contains a log-log plot of Pi as a function of 
Cbp. Two different regimes are clearly seen. First, in the 
low concentration region, the quantity Pi is nearly inde- 



pendent of Cbp. Second, in the high region, Pi decreases 
as a function of Cb p , being approximately proportional 
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FIG. 4. Schematic illustration of the interaction be- 
tween a reactive species and DNA molecules. We as- 
sume that a DSB is produced with a fixed probability q, pro- 
vided that at least one DNA molecule exists in the sphere of 
radius P rca + R K . 




FIG. 5. Theoretical prediction for the dependence of 
the threshold DNA base-pair concentration on the 
total length. This phase diagram in the Z-Cb P plane was 
obtained using equation JBJ. Region A: Pi is independent of 
Cbp. Region B: Pi is dependent on Cb p . 

The interaction of DNA molecules with ionizing radia- 
tion can be separated into two categories f§-fllj|. gdirect 
effects" (from energy deposited in the DNA and hydrated 
water molecules) and gindirect effects" (from free radi- 
cals produced by energy deposited in water molecules and 
other biomolecules located close to the DNA). In dilute 
solution, the ratio of the energy deposited in DNA and 
hydrated water molecules would be much smaller than 
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that deposited in the other molecules of the solution. In 
addition, Ito et al. has shown that the contribution of 
indirect effects to DSBs at Cb p = 100 /xM is approxi- 
mately two orders of magnitude larger than that of di- 
rect effects ■ Therefore it is reasonable to assume that 
indirect effects play a dominant role in the production of 
DSBs in our experiment. 

Now we propose a simple model to account for our ex- 
perimental data. In this model, it is assumed that each 
7-ray photon deposits energy along the radiation track 
mainly through Compton scattering (see FIGH]) [l_2j and 
produces reactive species [l3j]. Each reactive species has 
an effective diffusion length i? rea - According to a scaling 
argument, the radius of gyration R g can be regarded as 
defining the dimensions of a real DNA chain in a good 
solvent Given that at least one DNA molecule ex- 
ists in the sphere of radius i? rea + Rg , we assume that a 
DSB is produced with a fixed probability q (see FIG[4|). 
The sphere positioned at the center-of-mass of a reactive 
species means the gsphere of accessibility" for that reac- 
tive species : It is capable of attacking any DNA molecule 
located within in this sphere. Based on our framework, 
Pi is given by 



Pi = 



N Ica q{l - Q) 
NbpNnwA 



(1) 



where N rca is the number of reactive species produced 
per unit Gy, Q is the probability that no DNA molecule 
exists in a given sphere, iVb p is the number of base pairs 
for a single DNA molecule, and -/Vdna is the number of 
DNA molecules in the sample of total volume Vtot- The 
quantity Q is obtained by assuming that a single DNA 
moelcule uniformly distributes in the region, that has the 
radius R lcei + R g and the volume V, in the sample of total 
volume Vtot and that the positions of all DNA molecules 
are mutually independent variables. This yields the fol- 
lowing: 



Vtot 



Nona 



Next, because -/Vdna is given in terms of Cb p as 

^totCbp 



DNA 



A^bp ' 



for V <C Vtot, we obtain 



Then, combining equations (Q}, ([3]) and ((4]), we find 
iV rca g(l-e"^ Cbp ) 



Pi 



Vtot Cbp 



(2) 



(3) 



(4) 



(5) 



We have fitted our data with equation §5§ by treating 
N rea ,q/Vt t and V/N Dp as fitting parameters (see FIG. 



©). The best fit is obtained with N rca q/V to t = 1.31 /xM 
(lOOkbp-Gy)- 1 and V/N hp = 1.38 (juM) -1 . In addition, 
V/N hp is given in terms of i? rea and R g as 



V 
Abp" 



7r(i? rG 



bp 



(6) 



Then, substituting the best-fit value for V/ Nh P and the 



values N hp — 166 x 10 3 and R g = 1.5 /im lj|[l6j], we ob- 
tain i? rca ~ 1/im from equation ([6]). Finally, the effective 
lifetime of the reactive species, r, is estimated to be on 
the order of one millisecond from the relation i?^ oa ~ Dt, 
where D (~ 10~ 9 m 2 /s) is the diffusion coefficient fl7j . 
This value of r is reasonable, because the reactive species 
exist for approximately a millisecond (18(. 

The functional form of equation ([5]) suggests that the 
reshold cor 
ing relation: 



threshold concentration C£ p can be defined as the follow- 



C, 



bp 



iVbp 
V ' 



(7) 



The concentration CjJ p corresponds to the situation 
in which each region of volume V contains one DNA 
molecule and each reactive species attacks just one DNA 
molecule, on average. Therefore, the dependence in the 
high concentration region can be understood from the 
conjecture that the number of reactive species that can 
possibly attack a DNA molecule is essentially indepen- 
dent of Cbp, and thus increasing Cb p results in the re- 
duction of such reactive species per DNA molecule. 

Furthermore, we take into account the size of DNA in 
order to study its dependence on C^ p . From equations ([5]) 

and_0 and the relations N hp = f and R g = §A^ /5 Z 3 / 5 



C{J is obtained as follows: 



C, 



bp 



7r(i? r 



ca T 3 A k n ) 



(8) 



where rj (=3.4 A) is the distance between adjacent base 
pairs, Ak is the Kuhn length, and Z is the total length 
of a single DNA molecule. The quantity C^ p given by 
equation ([SJ is plotted as a function of Z in the case 
that i? rca = 1 /im and = 0.1 /xm in FIGJ5] This figure 
represents a phase diagram of Cb p and Z: In region A, Pi 
is independent of Cb p , while in region B, Pi is inversely 
proportional to Cb p . 

We now consider the process in which a single DSB is 
induced by reactive species. From the relations Pq (x I 
(see FIGE]) and I cx iVpho (the number of irradiating pho- 
tons), we obtain the relation Pq oc iV p ho- This assumes 
that DSB events produced by different photons are inde- 
pendent. Each gamma-photon produces clusters of reac- 
tive species along the radiation track 13, 111 ■ The diame- 
ters of these clusters are several nanometers, and they are 
widely separated in the case of radiation with low linear 
energy transfer 13, 19[ ■ Futhermore, as shown in FIGJ3J 
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Pi is inversely proportional to Cb p above ; i.e. the 
total number of DSBs is independent of Cb p . This im- 
plies that all reactive species attack DNA molecules and 
subsequently all of them disappear essentially. Therefore, 
this result indicates that the probability that a DSB will 
be produced is proportional to the number of reactive 
species that attack a single DNA molecule. With the 
scenario described to this point, we have two possibili- 
ties with regard to how a single DSB is produced: One 
possibility is that a single reactive species produces a 
single DSB with some probability, and the other is that 
a single cluster diffuses as an aggregate, attacks a sin- 
gle DNA molecule and produces a single DSB with some 
probability. However, the second possibility is seen to 
be unlikely, because the reactive species of each cluster 
diffuse almost homogeneously [l3|, [l^|. This argument 
thus leads to the conclusion that a single reactive species 
induces a single DSB with some probability (see FIG|4}. 
Our model, which is based on this scenario, yields results 
consistent with our experimental data(see FIG[3]). 

In conclusion, employing the single-molecule observa- 
tion method, we find that Pq is linearly dependent on I 
less than 200Gy, and that Pi remains nearly constant be- 
low a threshold value C£ and P\ is nearly inversely pro- 
portional to Cbp above C£ . Our theoretical model indi- 
cates that Pi , the number of DSBs per base pair per unit 
Gy, decrease as the DNA concentration increases above 
a certain threshold because the number of potentially at- 
tacking reactive species per DNA molecule also decreases. 
We conclude that it is necessary to consider not only the 
microscopic properties of the DNA molecule, but also 
the characteristics of giant DNA molecules as semiflex- 
ible polymers and the density characterizing the entire 
collection of such molecules for the systems in which they 
exist. This is a significant finding in our endeavor to un- 
derstand the physics of radiation-induced ionization of 
DNA molecules, and it also has important implications 
for biology and medicine. 



METHODS 

Experimental procedures and evaluation 
method. We used T4 phage DNA (166kbp, Nippon 
Gene) in a 20 mM phosphate buffer (pH 7.3) solvent. 
DNA was irradiated by 7-rays from a Co 60 source at 
a dose rate of 0.002-7 Gy/min. The fluorescent dye 
YO-YO-I (Abs/Em 491/509 nm, Molecular Probes, 
USA) was used for visualization of the DNA chains. 
After the irradiation, we added YO-YO-1 at 1/iM and 
2-mercaptoethanol (Nacalai Tesque, Japan) at 4% (v/v) 
at the final concentration in order to avoid additional 
oxidative damage. To fix the DNA molecules on a glass 
substrate, a glass base dish of 35 mm diameter (Iwaki, 
Japan) was cationically modified with a 0.05% (v/v) 
poly-L-lysine solution (Sigma- Aldrich, USA), and a 



sample solution (~5 /A) was adsorbed onto the modified 
glass base dish. To elongate the DNA molecules, the 
droplet of sample solution was covered with a glass 
slip (18 mmxl8 mm). We applied a weak shear to 
the DNA molecules by sliding this glass slip. This was 
done very gently to avoid additional fragmentations 
by mechanical stress. Fluorescence images of DNA 
molecules were observed using an Axiovert 135 TV (Carl 
Zeiss, Germany) fluorescence microscope equipped with 
a 100 x oil-immersed objective lens. The lengths of the 
elongated DNA molecules were measured using cosmos 
image- analysis software (Library, Japan). We measured 
the contour lengths of approximately 100 molecules 
in each sample solution and thereby obtained L(I). 
Giant DNA molecules exceeding lOOkbp are occasionally 
inadvertently cut into fragments during the sample 
preparation, which includes pipetting and mixing |20j. 
For this reason, the value of L(0) is sli ghtl y shorter than 
the natural length of T4 DNA, 57 //m [El, H3|. 

The quantity Po (the number of DSBs per base pair) is 
calculated from the difference between the average molec- 
ular lengths of irradiated and non-irradiated samples that 
are otherwise identically prepared: 



_ L(0)/L(I) 1 

p °- — W) — 



(9) 



Here, 77 (=3.4 A) is the distance between adjacent base 
pairs. 
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